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ABSTRACT 

We present a statistical analysis of simultaneous optical and X-ray light curves, spanning 
600 ks, for 814 pre-main-sequence (PMS) stars in the Orion Nebula Cluster. The aim of this 
study is to establish the relationship, if any, between the sites of optical and X-ray variability, 
and thereby to elucidate the origins of X-ray production in PMS stars. In a previous paper we 
showed that optical and X-ray variability in PMS stars are very rarely time-correlated. Here, 
using time-averaged variability indicators to examine the joint occurrences of optical and X-ray 
variability, we confirm that the two forms of variability are not directly causally related. However, 
a strong and highly statistically significant correlation is found between optical variability and 
X-ray luminosity. As this correlation is found to be independent of accretion activity, we argue 
that X-ray production in PMS stars must instead be intimately connected with the presence 
and strength of optically variable, magnetically active surface regions (i.e. spots) on these stars. 
Moreover, because X-ray variability and optical variability are rarely time-correlated, we conclude 
that the sites of X-ray production are not exclusively co-spatial with these regions. We argue 
that solar-analog coronae, heated by topologically complex fields, can explain these findings. 

Subject headings: open clusters and associations: individual (Orion Nebula Cluster) — stars: 
flare — stars: magnetic fields — stars: spots — stars: pre-main-sequence — X-rays: stars 



Introduction 



While it is now well established that low- mass pre-main-sequence (PMS) star s typically have X-ra y 



Preibisch ct al 



2005). 



luminosities 2-4 orders of magnitude higher than that of the present-day Sun (e.g. 
fundamental questions remain as to how these X-rays are actually produced. Observations of X-ray flares 
in these stars indicate solar-analog magnetic reconnection events, but the absence of an obvious main- 
sequence-like rotation-activity relation in PMS s tars has so far confounded attempts to ascribe PMS X-ray 



production directly to d ynamo-generated fields (|Flaccomio et al.ll2003t iFeigelson et al.ll2003t IStassun et al 



20041: iRebull et al.ll2006h . Accretion has also been implicate d as a possible driver of PMS X-ray produc- 
tion, but no consen sus has emerged on this point as yet (e.g. Kastner et al. 20021 : Stelzer fc Schmitt 2004 : 
Stassun et al.ll2004h. Indeed, there is evidence that accretion m ay act to inhibit X-ray production instead 
of enhance it (e.g. Flaccomio et al. 2003 : Preibisch et al. 2005h . These issues are reviewed more fully in 
Feigelson et all ( 200 1 
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X-rays from PMS stars are interesting not only because of the mystery surrounding their origins. If 
X-rays trace magnetic activity in young stars as they do in main sequence stars, then X-rays can provide 
important clues to a variety of star-formation questions involving stellar magnetic fields — the evolution of 
stellar angular momentum via magnetized winds and magnetic coupling to circumstellar disks, for example. 
X-rays may also be central players in the phys ics of molecular clounds and protoplaneta ry disks, driving the 
photoionization of circumstellar material (e.g. Glassgold. Feigelson. fc Montmerlel 2000f ). Thus, understand- 



ing the physics of X-ra y production in PMS stars may ultimately prove central to our understanding of star- 
and planet-formation ( Feigelson et al. 200(jh . 



To make progress in understanding how PMS X-rays are made, it would be valuable to have constraints 
on where the X-rays are coming from. A key limitation in this regard has been the lack of simultaneous, 
independent measures of X-ray production and of other forms of activity, such as optical variability. So 
motivated, we have undertaken an extensive variability study of the Orion Nebula Cluster (ONC) combining 
a nearly continuous, 13-d Chandra observation with simultaneous, multi-wavelength, time-series photometry 
in the visible. The resulting database of simultaneous X-ray and optical light curves for some 800 members of 
the ONC represents, by a factor of hundred s, the largest attemp t to date to study the relationship betweeen 



X-ray and optical variability in PMS stars (jStassun et al.ll2006l hereinafter Paper I) 



In Paper I, we used this database to search for time-correlated variability in our simultaneous optical 
and X-ray light curves with the goal of establishing the frequency with which X-ray and optical variations 
are co-temporal, and to thereby constrain the extent to which the sources of X-ray and optical variability 
may be co-spatial or otherwise causally linked. For example, if X-rays are produced near accretion shocks 
on classical T Tauri stars (CTTSs; actively accreting PMS stars), then one might expect that changes in 
the strength of these shocks will induce changes in the strength of the X-ray emission, and that X-ray and 
optical variability might therefore be correlated in time. Similarly, if X-rays originate in coronal structures 
associated with magnetically active surface regions on these stars (i.e. dark spots), we might then expect 
that optical and X-ray variability would be anti-correlated in time. However, we found that such behavior 
almost never occurs. Indeed, except for a single possible case of "white light" flare emission, simultaneous 
X-ray/optical variations were not clearly seen in any of the 800+ stars observed. The upper limit on the 
incidence of time-correlated optical/X-ray variability in our sample is at most ~ 5%; that is, at least ~ 95% 
of the stars in our sample exhibit no evidence for time-correlated optical and X-ray variability (Paper I). The 
implication is that the sites of optical and X-ray variability in PMS stars are not — in the vast majority of 
cases — instantaneously one and the same. Sim ilar findings have now also been reported from a simultaneous 



optical/X-ray monitoring survey of Taurus bv lAudard et al.l |2007) 



Whereas these findings rule out the simple interpretation of the X-ray and optical variability both 
dominated by one site of emission, in this companion paper we examine the possibility of a more general 
relationship between optical and X-ray variability by investigating time- averaged optical and X-ray vari- 
ability indicators (e.g. variability amplitudes) and possible correlations therein. After reviewing the salient 
contents of our simultaneous optical/X-ray variability database, and the time-averaged variability indicators 
that we use in our analysis ($21), m fj3] we perform a multivariate correlation analysis of these indicators 
with one another and with other stellar properties related to energy production (accretion signatures, bolo- 
metric luminosity, X-ray luminosity). We find that time-averaged optical and X-ray variability are not well 
correlated with one another, implying that they are not causally related. Importantly, however, we do find 
optical variability to be strongly correlated with X-ray luminosity, and that this correlation exists regardless 
of whether the stars in question are actively accreting (CTTSs) or non-accreting (weak-lined T Tauri stars; 
WTTSs). 
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We discuss in S}4]the implications of these results. We argue that accretion is not in general a significant 
contributor to the production of X-rays in CTTSs, and that PMS X-ray production overall is more intimately 
related to the presence of magnetically active surface regions (i.e. spots), as traced by optical variability, in 
both CTTSs and WTTSs. At the same time, the lack of time-correlated optical and X-ray variability in the 
overwhelming majority of our sample (Paper I) indicates that the X-rays are not produced in magnetically 
active surface regions alone. We summarize our conclusions in $5] 



2. The COUP synoptic X-ray and optical database 



A complete description of our database of simultaneous X-ray and optical light curves is provided in 
Paper I. Briefly, the database comprises 814 PMS members of the ONC within a 17'xl7' region centered on 
the Trapezium observed by ACIS onboard Chandra for 13.2 d. These stars were simultaneously observed at 
BVRI wavelengths by small telescopes on the ground during the last 7 d of the Chandra observation, with 
a cadence of ~ 1 hr _1 . 



For each star, the database also includes the star's X-ray luminosity, Lx ( Getman et al. 20051 ). and 
a compilation of other physical parameters fro m photometric an d spectroscopic catalogs already in the 
literature, including bolo metric luminosity, Lboi ( Hillenbrand! 1997 ). and the equivalent width (EW) of the 
Ca I I infrared triplet from Hillenbrand et al. ( 1 99q). a measure of mass accretion rate. As in previous studies 



(e.g. 



Flaccomio et al. 200 



3; 



Stassun et al 



2004 IPreibisch et aDl2005l . Paper I), we take the "accretors" to 
be those stars with Ca II strongly in emission, i.e., EW < — 1 A, and "non-accretors" to be those with Ca II 
clearly in absorption, i.e., EW > 1 A. Of the 493 stars in our study sample with Ca II measurements, 151 
stars are accretors and 145 stars are non-accretors, as defined here (the remainder show indeterminate values 
close to Afl 

The 814 stars in our database span a very large range of apparent magnitudes and colors, with 7.3 < 
/ < 18.4 and —0.1 < (V — I) < 5.7, and a correspondi ngly large range of spectral types and extinctions, B3 
< SpTy < M7 and 0.0 < Ay < 10.8 ( Hil lenbrandlll997n . This sample is thus representative of the underlying 



ONC population; only the very massive, heavily obscured, and sub-stellar populations are excluded. 

These stars exhibit a large variety of behaviors in their optical and X-ray variability. Figs. 1-6 in Paper I 
show the optical and X-ray light curves for a selection of objects as a representative visual summary of the 
range and types of variability observed. Below we describe the statistical indicators that we use to quantify 
the time-averaged X-ray and optical variability. 



1 As a caveat, we note that EW[Ca II] is an imperfect indicator of accretion. As shown by Isicilia- Ag uilar et alj ll2005h . some 
PMS stars with clear accretion signatures in Ha do not manifest this accretion activity in Ca II. It is therefore possible or 
even likely that some low-level accretors will go undetected in Ca II. The criterion that we adopt for identifying the presence 
of accretion is thus in practice a conservative one; stars with EW[Ca II] < -1 A and EW[Ca II] > 1 A are securely identified 
as accretors and non-accretors, while some stars with indeterminate values close to A may be weakly accreting objects that 
will go unidentified as such. 
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2.1. Time-averaged variability indicators 

2.1.1. X-ray variability 

The COUP database provides two measu res of X-ray variability in the detected X-ray sources. These 



are discussed in detail bv lGetman et al.l (|2005f ) 



The first measure is a standard one-sided Kolmogorov-Smirnoff (KS) test, which tests the null hypothesis 
that the photon arrival times are distributed uniformly in time. The KS test thus gives the probability that 
the source's X-ray flux is consistent with being non- variable (i.e. constant). We take sources with KS test 
probabilities of less than 0.001 as being definitively X-ray variable, the expected number of false positives in 
our sample then being less than 1. Of the 814 sources in our study sample, 556 (68%) are X-ray variable at 
greater than 99.9% confidence by this definition. 

The second measure of time-averaged X-ray variability is a Bayesian Block (BB) analysis ( Scarglel 



19981 ). which segments each light curve into the maximum number of time blocks such that the differences 
in the mean flux levels of the blocks are statistically significant. The BB analysis thus yields a robust 
measure of the maximum and minimum flux levels in the source's X-ray light curve (BB max and BB m j n , 
respectively). In our analysis, we use the ratio BB max /BB m i n to quantify the magnitude of variability in the 
X-ray light curves of stars in our sample. The X-ray variables in our database have a range of amplitudes 
0.12 < log(BB max /BB m i n ) < 2.89 (see Fig.Q]). The KS and BB measures of time-averaged X-ray variability 
for the stars in our study sample are reported in Table [TJ 



2.1.2. Optical variability 

To char acterize the tim e-averaged optical variability of the stars in our sample we use the J statistic 
described by Stetson ( 1996h . Similar to the standard x 2 statistic, the J statistic measures the residuals of 



the data about the weighted mean level of the light curve, and compares those residuals to the measurement 
uncertainties. However, unlike a standard \ 2 analysis, the Stetson index simultaneously considers the mea- 
surements from light curves in all available photometric bands of a given source. This is done by pair-wise 
matching the data points in time and comparing the signs of the matched residuals. 

j= ^ =1 w t sgn(P l )^\P~~\ 



where p is the number of pairs of observations of the star, Pj = <5jYi)^fcr»i is the product of t he normalized 
residuals of the paired measurements, and wi is the statistical weight assigned to each pair (see Stetson 19961 
for details on weights). The normalized residual of measurement i in passband j is given by 



5 



n - 1 



where n is the number of measurements used to determine the mean magnitude m, and Oi is the photometric 
uncertainty in measurement m^. 

The J statistic is more robust than \ 2 against outliers. As an example, consider a non-variable star 
with simultaneous light curves in three passbands, and consider the case where a cosmic ray hit appears 



- 5 - 



as a discrepantly high data point in one of these three light curves. This one outlier will increase the 
value of the \ 2 statistic for that light curve, making it appear variable in one passband but not in the 
others. I n con trast, this outlier will tend to decrease the value of the J statistic. Following the procedure of 
Stetson ( 1996f ) to make the J statistic even more robust, we compute the J statistic iteratively, adjusting the 



weights of individual data points and re-computing the weighted mean of each light curve between successive 
iterations, thereby minimizing the influence of strong outliers on the J statistic. 

To calculate J in practice, the light curves must be binned in time, which requires a choice of binning 
timescale. This binning is not in the sense of combining adjacent data points; rather, time bins are defined to 
determine which data points from the multiple light curves are associated with one another in the pair-wise 
matching discussed above. Given the typical sampling and cadence of the optical light curves (see Paper I), 
we adopt time bins of 1 hr. 

Values of J for the stars in our study sample are reported in Table [TJ For each star we also report in 
Table [1] the measures of time-averaged X-ray variability and of accretion discussed above. The distribution 
of J for our study sample is shown in Fig. [TJ From this distribution, we can empirically choose a threshold 
of J, above which we may be confident that the optical light curves manifest true variability. True variable 
stars will show positive values of J, the magnitude of which is related to the amplitude of the variability. 
In comparison, non-variable stars will produce a narrow distribution of J centered about zero, with equal 
probability for being positive or negative. This narrow distribution of values about zero represents the 
probability distribution of J arising from small, chance correlations in the optical light curves, and we may 
thus use the most negative value of J as a conservative estimate of the most positive value of J that will 
result from such chance correlations. 

The distribution of J in our sample (Fig. [TJ) shows a maximal negative value of J = —1.50. This value 
occurs once in our sample of 814 stars, and we may therefore expect that a value of J > 1.50 will occur by 
chance with a probability of 1/814 rj 0.001, which is the same as the conservative criterion we have adopted 
for identifying variability in the X-ray light curves. 

With this threshold, 358/814 (44%) of the stars are optically variable, with a range 1.5 < J < 45.5. 
The optical and X-ray light curves for a selection of sources with J spanning this range are shown in Figs. 
1-6 in Paper I. In general, larger values of J imply greater variability. The threshold of J > 1.50 that 
we adopt for identifying truly variable objects corresponds roughly to peak-to-peak variability of ~ 0.05 
mag (typical photometric errors are 0.01-0.02 mag; Paper I), though for any particular source the peak-to- 
peak variability is different at the different wavelengths (generally larger at shorter wavelengths due to the 
increased temperature contrast). Thus, sources identified as "non- variable" in fact represent upper limits of 
peak-to-peak variability less than ~ 5%. At the other extreme, objects with the largest J values show peak- 
to-peak variability in the V band of > 2 mag. With this range of variability amplitudes, the time-averaged 
optical var iability of our samp le is typical of what has been observed before in optical variability studies of 



TTSs fe.g. lHerbst et al.lll994h . 



3. Results: Statistical analysis of time-averaged variability indicators 

Here we report the results of a statistical analysis of these variability indicators, with the goal of 
determining their relationships with one another and with other stellar parameters. We present three levels 
of analysis in order of increasing sophistication ( ^3.1H3.3p with a concise summary in £13.41 All statistical 
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analyses reported here were performed with the R software^ package 



3.1. Relation of optical and X-ray variability to accretion 

One of the aims of our study is to determine whether accretion may be a significant contributor to the 
production of X-rays in PMS stars. We have already seen in Paper I that there is no evidence in the form 
of time-correlated X-ray and optical variability for such a connection, implying that any causal relationship 
between accretion and X-ray production must be indirect. Here we consider whether the presence of active 
accretion manifests itself in the time-averaged variability of our sample. 

Applying a Student's t test to the distributions of X-ray variability (BB max /BB m i n ) for the accretors and 
non-accretors, we find a probability of 38% that these two groups have similar mean values of BB max /BB m ; n . 
That is, there is not a statistically significant difference in the level of X-ray variability among the accretors 
as compared to the non-accretors. 

In contrast, accretion does manifest itself in the optical variability of these stars. A t test applied to the 
same accretor and non-accretor samples confirms, at the 99% confidence level, that accretors exhibit higher 
average levels of optical variability (J) than the non-accretors. This is as expected, since optical variability 
in CTTSs is believed to be driven in part by accretion. This result is also reassuring, as it demonstrates 
that the effects of accretion are discernible in the gross statistical behavior of our study sample, even though 
multiple processes may be contributing at once to the observed variability, and even though the EW[C a II] 
measurements that we use to identify the accretors were obtained years earlier ( Hillenbrand et al. 1998t h 



3.2. Joint occurrences of optical and X-ray variability 

Having found that accretion manifests itself in the level of optical variability of our sample, but not 
in the level of their X-ray variability, we next compare the incidences of optical and X-ray variability in 
these stars. We find statistically significant optical variability among 358/814 (44%) of the stars (Table 
[2]) and statistically significant X-ray variability among 556/814 (68%) of the stars ( *12.1.ip . A Pearson % 2 
contingency-table test shows that the probability that optical and X-ray variability occur with the same 
frequency is less than 10~ 6 . In other words, a statistically significant fraction of stars that exhibit X-ray 
variability do not also exhibit optical variability in our data. 

Conversely, we consider whether optically variable stars always exhibit X-ray variability. Of the 358 
stars in our sample that we have identified as optically variable, 278 (78%) of them also exhibit X-ray 
variability (see Table [J). The likelihood that this fraction is consistent with 100% is again less than 10~ 6 . 
We can therefore state that, to the sensitivity available, the cause of optical variability cannot at the same 
time be the cause of X-ray variability for a significant fraction of the 814 observed stars. The two forms of 
variability simply do not always occur together during the 7 days of overlap between the optical and Chandra 
observations. 

Perhaps the occurrences of optical and X-ray variability are more strongly related among those stars 
that are accreting, as one might expect if active accretion drives both forms of variability. To test this, 
we again segregate the stars in our sample into accretors and non-accretors, as defined earlier. Among the 



2 R is a freely available language and environment for statistical computing and graphics. See http://www.r-project.org 
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accretors, we find simultaneous occurrences of optical and X-ray variability for 60 stars (40%), whereas for 
the non-accretors we find simultaneous occurrences of optical and X-ray variability for 64 stars (44%) (see 
Table [5]). This small difference is not statistically significant (52% probability that they are not different) 
according to a Pearson \ 2 contingency-table test. 

While we find clear evidence that optical and X-ray variability are not always linked to one another, and 
that the presence of active accretion does not increase the likelihood for an association between the two, this 
does not necessarily imply that the occurrences of optical and X-ray variability are completely independent 
of one another. In fact, we find that stars showing one form of variability are more likely to also show 
variability in the other. For example, whereas optical variability is found in 44% (358/814) of the overall 
study sample, this fraction increases to 50% among the stars in our sample that are X-ray variable (278/556), 
and 69% (178/258) of the X-ray non- variable stars are also optically non- variable (see Table [2|). A Pearson 
X 2 contingency-table test shows that the likelihood of the two measures of variability being independent 
of one another is less than 10~ 6 . Thus, optical and X-ray variability are, in fact, related. However, this 
relationship is not directly causal, as we now show. 



3.3. Strength of optical and X-ray variability and their relation to X-ray luminosity 

To explore the possible link between optical and X-ray variability further, we next examine how the 
strength of one may be related to the other. In Fig. [21 we plot the strength of optical variability, as measured 
by the J statistic f ^2.1.2[) . versus the strength of X-ray variability, as measured by the BB analysis f £|2.1.ip . 
A Kendall's r test shows that the two measures are correlated at a high level of statistical significance, 
with the probability that they are not correlated being 3 x 10~ 6 . The same test conducted separately 
among the smaller sub-samples of accretors and non-accretors shows a correlation in the same sense, albeit 
at marginal levels of significance (Kendall's r probabilities of 0.08 and 0.009 for the accretors and non- 
accretors, respectively). Thus, we find that not only does the presence of one form of variability correlate 
with the other (see above), but that optical and X-ray variability are furthermore correlated in strength. 

But do these statistical correlations between optical variability and X-ray variability imply a causal 
relationship between the two? The data suggest that this is not the case, and that the statistical relationship 
between optical and X-ray variability is in fact due to a mutual correlation — the first physical, the second 
observational — with X-ray luminosity. 



X-ray luminosity (Lx) was found bv IStassun et al.l (|2004f ) to be strongly correlated with optical vari 



ability in their analysis of all previous Chandra observations of PMS stars in the ONC. Specifically, they 
found that the ONC stars with the highest Lx also exhibited the highest levels of optical va riability. This 



correl ation is confirmed here, with a sample that is five times larger than that studied by IStassun et al 



(|2004l ). The correlation between J and Lx in our study sample (Fig. [3^) is highly statistically significant: 
The probability that J and Lx are uncorrelated, according to a Kendall's r test, is less than 10~ 6 . Impor- 
tantly, we furthermore find that this correlation persists, at a less strong but still highly significant level, 
when we separately consider accretors and non-accretors. 

At the same time, our sample exhibits a similar correlation between Lx and X-ray variability (Fig. [3Jd), 
with similarly high statistical significance as that between Lx and optical variability. In this case, however, 
the correlation appears to stem from an observational bias against detection of X-ray variability in stars with 
low Lx- Indeed, as shown in Fig. [3j3, all of the stars for which the BB analysis finds no X-ray variability 
(BB max /BB min = 1) have L x < 10 30 erg s" 1 . 
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That optical and X-ray variability are both so strongly correlated with Lx suggests that this mutual 
correlation may be responsible for the correlation of optical and X-ray variability with another. To bet- 
ter understand these mutually correlated variables, we perform a multiple linear least-squares regression 
analysis on the strength of optical variability (J) simultaneously against the strength of X-ray variability 
(BB max /BB m i n ) and Lx- The regression coefficient for Lx is found to be statistically significant (null hy- 
pothesis probability of 5 x 10 -6 ), but the coefficient for BB max /BB m i n is not (null hypothesis probability of 
0.34). Similarly regressing BB max /BB m i n simultaneously against J and Lx, we find that the Lx coefficient 
is statistically significant (null hypothesis probability of less than 10 -6 ) but that the J coefficient is not (null 
hypothesis probability of 0.34). The regressions of J and BB max /BB m i n vs. Lx are represented by dashed 
lines in Fig. [3] 

Thus, the relationship between optical and X-ray variability in our sample is, in fact, fully explained 
by the mutual correlations of J and BB max /BB m i n with Lx- Controlling for this mutual correlation with 
Lx, we find no statistically significant residual correlation of optical and X-ray variability with one another. 
Optical and X-ray variability are evidently not causally related. 

Even though the correlation of optical variability with Lx has a plausible phys ical explanation (stars 



that are more heavily spotted are also stronger X-ray emitters; IStassun et al.l 120041 ). our finding that the 
correlation of X-ray variability with Lx is due to observational bias (Fig. [HJd) raises the suspicion that the 
correlation of optical variability with Lx may be similarly spurious. For e xample, Lx is known to be very 



strongly correlated with bolometric luminosity (Lboi) m the COUP sample (jPreibisch et al.ll2005| ). If optical 
variability is in turn correlated with Lboi in a manner similar to the correlation of X-ray variability with Lx , 
then this may explain the correlation of optical variability with Lx as being non-physical. 

To examine this possibility quantitatively, we again perform a multiple regression analysis, as above, 
but this time on J against Lboi and Lx- First of all, we find that optical variability is indeed correlated with 
Lboi, albeit weakly (null-hypothesis probability of 0.02), when we regress J against Lboi alone. However, 
something unexpected occurs when we include Lx as a regressor. The regression coefficient for Lboi becomes 
statistically insignificant (null- hypothesis probability of 0.76), whereas the coefficient for Lx is significant, 
and with a higher degree of significance (null-hypothesis probability of O.OO^jfl) than was the correlation 
of J with Lboi alone. In other words, controlling for the correlation of J with Lboi does not eliminate 
the correlation of J with Lx', on the contrary, controlling for the correlation of J with Lx eliminates the 
correlation of J with Lboi- Lx is evidently the more fundamental correlate. 

This somewhat surprising result reinforces the conclusion that optical variability is intrinsically related 
to X-ray luminosity in our study sample. Having considered in detail the mutual correlations of J with 
BB max /BB m j n , Lboi, and Lx, only the correlation between J and Lx remains. The correlation between 
optical variability and Lx is more significant even than the intuitive correlation of optical variability with 
accretion f ij3.f p . Indeed, after controlling for the influence of all other secondary variables, we find that the 
strongest correlate of optical variability in our study sample is X-ray luminosity. 



3 Note that the statistical significance of the correlation between J and Lx is lower here than in the multiple regression 
analysis of J on Lx and BB max /BB m ; n . This is because the sample under consideration here is smaller (545 stars vs. 814 
stars), since Lboi measurements are not available for all stars in the COUP database. 
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3.4. Summary 

We find that optical and X-ray variability in our study sample are not causally related. This conclusion 
follows from an examination of the joint occurrences of time-averaged optical and X-ray variability, as well 
as their mutual correlations with one another and with other variables. We find definitive evidence that, on 
timescales of 7 d or less, optical variability does not always occur in the presence of X-ray variability, nor 
vice-versa. Even when the two forms of variability are present, non-accreting stars are just as likely to show 
both forms of variability as are stars that exhibit active accretion; that is, among stars exhibiting both optical 
and X-ray variability, accretion is evidently not the cause. At a more detailed level, we have performed a 
multiple regression analysis, controlling for the mutual correlations of optical and X-ray variability, both 
with one another and with bolometric and X-ray luminosity. This analysis reveals that, while statistically 
significant correlations exist among all of these measures, the fundamental correlation in this study is that 
between optical variability and X-ray luminosity. This correlation is furthermore independent of accretion, 
indicating again that accretion is not the cause. 



Discussion 



While magnetic activity, and therefore X-ray produc tion, in late-type main-se quence stars is relatively 
well understood in terms of a rotation-driven dynamo (e.g. lSchriiver fc Zwaa whether stellar rotation 

is also at the heart of PMS X-ray production is a contested iss ue. In contrast to the strong "rotation-activity 
relationship" that is clear ly present on the main sequence (|Pallavicini et al.l Il981t Ijeffriesl 11999c iRandich 
2000 ; Pizzolato et al. 2003h . a series of studies based on Chandra observations of a large sample of PMS stars 
in the Orion Nebula Clus ter (ONC) failed to find any direct correlation between the s t ars' rotation rates an d 
their X-ray luminosities ( Flaccomio. Micela. fc Sciortino 2003bl ; Feigelson et al. 2003 ; Preibisch et al.| [2005). 



However. IStassun et al.l (|2004r ). in an analysis of all archival Chandra observations of PMS stars in the 
ONC, found indirect evidence hinting that a main-sequence- type rotation- activity relationship may in fact be 
driving PMS X-ray production after all. They found a statistically significant, positive correlation between 
stellar rotation period and Lx, suggesting that these stars are on the "super-saturated" portion of the 
rotation-activity relation, consistent with their Rossby numbers. They furthermore found that the sample of 
ONC stars with known rotation periods is significantly biased to high Lx] stars whose rotation periods are 
unknown, and which have lower L x , may thus repr e sent t he unseen "linear" portion of the rotation-activity 
relation. Preibisch et al. ( 2005h and Rebull et al. (2006) have since confirmed these findings in different 
Orion samples. 

The results of the present study lend further support for an interpretation in which accretion is not a 
significant source of X-ray production in PMS stars. Accretion, as traced by emission in Ca II, does not 
manifest itself in the X-ray variability of our study sample ( §3.ip . And while accretion is manifested in the 
optical variability of these stars as expected, a multiple regression analysis finds that X-ray variability is not 
correlated with this optical variability f ^3.3|) . Moreover, we have found very few cases in our study sample 
in which the optical and X-ray variability are correlated in time (Paper I). 

That we find so little evidence for accretion as a dominant source of X-rays in young stars suggests 
that we look to other possible X-ray production mechanisms. Indeed, there are now strong indications 
that X-ray production is intimately related in some way to the presence of magnetically active surface 
regions on these stars. The evidence for the existence of magnetically active surface regions on PMS stars is 
abundant. Several important studies of Zeeman broadening in the photospheric spectral lines of PMS stars 
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Guenther et al. 19991: 


Johns-Krull, Valenti. & Saar 2004) 


. Rice & Strassmeier 


1996: Johns-Krull & Hatzesl 1997). 



of the existence of strong surface fields. Moreover, the X-ray flares obs erved in the COUP s tudy can only 
be reasonably explained as reconnection events in strong stellar fields (jFeigelson et al.l 12006). In addition, 



numerous photometric studies of PMS stars have relied upon the presence of surface f eatures to measure 



these stars' rotation periods (e.g. IStassun et al.l Il999t iRebull! l200lh iHerbst et al.l I2002T ) . The stability of 



these features over timescales of up to several years strongly suggests that these features are in many cases 
"starspots" presumably analogous to dark spots on the Sun. 



Stassun et al.l (|2004f ) found that PMS stars in the ONC exhibit a highly statistically significant corre- 
lation between the strength of their optical variability (taken to be a measure of the stars' "spottedness" ) 
and their X-ray luminosities, in the sense that the least X-ray luminous stars are those with the smallest 
amount of spottedness. This finding led those authors to infer a causal connection between X-ray emission 
and the presence of starspots, leading to the conclusion that solar-analog surface magnetic activity driving 
hot coronae is the primary source of X-ray production in PMS stars. However, the possibility remained that 
the correlation between optical variability and X-ray luminosity was driven by accretion. 

The results of the present study firmly corroborate the correlation between optical variability and X- 
ray luminosity (i.e. Fig. [3]) — indeed, optical variability is the strongest correlate of X-ray luminosity in the 
COUP study sample (SJ3T3J and furthermore show that this correlation exists regardless of whether the 
stars under consideration are accreting or non-accreting ( fl3.ip . This finding contradicts models in which 
the observed X-ray emission is produced primarily by accretion, and instead bolsters support for the role of 
non-accretion-related, optically active magnetic regions (i.e. spots) in the production of PMS X-rays. 

A small number of the stars in our sample evi nce periodic X-ray va riability with a periodicity that is 
very similar to the observed periods in the optical ( Flaccomio et al. 2005h . In these cases the sites of X-ray 
emission may be predominantly co-spatial with magnetic surface spots. However, the association of X-ray 
production with magnetic spots does not manifest itself in a time-correlated way for most stars in our sample 
(Paper I), indicating that the X-ray emitting coronae in PMS stars generally have spatial structures that 
little reflect the spatial distribution of magnetic footpoints in the underlying photosphere. 

These conclusio ns are consistent with recent studies of TTS magnetic field structures. For example, 
Jardine et al. ( 2006h find that a variety of field topologies are required to explain the observed scatter of 
emission measures of COUP stars, ranging from simple dipole fields to much more complex configurations. 
As discussed by those authors, the more complex field configurations are characterized by coronal gas that is 
confined to compact loops covering a large fraction of the stellar surface. Coronal emission arising from such 
a topology will have a relatively low X-ray luminosity by virtue of its compact size, and the footpoints of 
the multi-polar field will produc e correspondingly small spots that are more-or-less uniformly distributed on 
the stellar surface (cf. Fig. 8 in I Jardine et al.ll2006l ), resulting in a relatively low- amplitude spot signal in the 
optical. The reverse will be true for a star whose coronal emission arises in a more extended, more dipolar 
field, rooted in a small number of footpoints (spots) that are more asymmetrically distributed on the stellar 
surface; such a configuration will produce a high X-ray luminosity together with strong optical variability. 
Thus, given a large ensemble of PMS stars with a range of field configurations, a strong correlation between 
X-ray luminosity and optical "spottedness," as we have found here, may naturally follow. 



5. Conclusions 



We have presented a combined analysis of time-averaged optical and X-ray variability for more than 800 
pre-main-sequence (PMS) stars in the Orion Nebula Cluster (ONC) with the aim of elucidating the origins 
of X-ray production in PMS stars. This study complements and extends the analysis of time-correlated 



variability in this same sample reported in lStassun et al.l ( 20061 ) 



Combining the results of these two studies, we find very little evidence to suggest a direct physical link 
between the sources of optical and X-ray variability in these stars. Whether we consider time-correlated 
behavior between the individual stars' simultaneous optical and X-ray light curves, whether we consider the 
optical and X-ray variability properties of our study sample as a statistical ensemble, whether we segregate 
the stars according to their accretion properties, we arrive at the same conclusion: the sites of optical and 
X-ray variability in PMS stars are not — for the vast majority of stars — instantaneously one and the same. 
An interpretation in which accretion is a dominant source of X-ray production in PMS stars is not supported 
by the data. 

Instead, the data provide strong support for an interpretation in which X-ray production in PMS stars is 
intimately connected with the presence and strength of optically variable, magnetic regions on their surfaces 
(i.e., spots). Indeed, after controlling for the influence of all other secondary variables, including accretion, 
we find that the strongest correlate of X-ray luminosity in our sample is spottedness (as measured by optical 
variability, after controlling for accretion). Our interpretation is that solar- analog magnetic activity is the 
primary, common driving mechanism of both optical variability and X-ray production in these stars. The 
fact that this relationship between spottedness and X-ray production does not, in the vast majority of cases, 
manifest itself in a time-correlated way (Paper I) suggests that the optical sp ots represent the foo tpoints of 
complex magnetic topologies that heat solar-analog coronae, as suggested by Jardine et al. ( 20061 ). 
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Tablc 1. Optical and X-ray Variability in Study Sample 
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Tablc 1 — Continued 
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Table 1 — Continued 
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a Optical identification from Gctman ct al. (2005). 
b Valuc of J statistic; sec 92.1.21 

C KS test probability of X-ray variabili ty; see ^^JLjJ 
Value s below —4.00 have been truncated jGetman et all 
12005ft . 

d Ratio of maximum/minimum X-ray flux from BB 
analysis; sec jj2. 1. 1| 

c Equivalent width in A of the Ca II IR triplet lines. 
Negat ive values indicate emission. See [Hillenbran d et al.| 
d 1998ft . 



Note. — This table is available in its entirety in the 
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electronic version of the Journal. A portion is shown here 
for guidance regarding its form and content. 
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Tabic 2. Contingency table of optical and X-ray variability in the study sample 



Optical X-ray variability 



variability Yes No 



Entire study sample (N — 814) 

Yes 278 80 

No 278 178 

Actively accreting stars (N — 151) 

Yes 60 25 

No 39 27 

Non-accreting stars (N = 145) 

Yes 64 8 

No 59 14 
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Fig. 1. — (a) Histogram of the J statistic for all stars in the study sample. Chance correlations in the light 
curves produce a distribution of positive and negative J values centered on zero. The most negative value 
of J is —1.50; thus J > 1.50 (indicated by the vertical line) serves as a conservative threshold for identifying 
stars that are likely true variables. Stars with various values of J above this threshold are indicated; their 
light curves can be found in Paper I. Inset panel compares the distribution of J separately for accretors and 
non-accretors as defined in ^ (b) Same as above but for log BB max /BB m i n . The first bin is dominated by 
stars with log BB max /BB m j n = 0, which are non-variable in X-rays. 
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Fig. 2. — Strength of optical variability as measured by the J statistic is plotted vs. strength of X-ray 
variability as measured by the BB analysis. Stars with negative J have been excluded from the plot so that 
it may be displayed on a logarithmic scale (see Fig. QJl for explanation of negative J values). The dashed 
line shows a linear regression fit to all of the data. A statistically significant correlation in the sense of the 
line is confirmed by a Kendall's r test. This correlation is found to be the result of the mutual correlation 
of J and BB max /BB min with X-ray luminosity; see text. 
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Fig. 3. — (a) Strength of optical variability as measured by the J statistic is plotted vs. X-ray luminosity 
{Lx )• (b) Strength of X-ray variability as measured by the BB analysis is plotted vs. Lx- Symbols and lines 
are as in Fig. O In both panels, a statistically significant correlation in the sense of the line is confirmed 
by a Kendall's r test. This correlation is in the same sense when the accretors or non-accretors are tested 
separately. Only the correlation in (a) is likely to be physical; see text. 



